4*0*11 QfrUll 

rtmihioUJIj rtmihiU 
griUrf 

Solving Homogenous and Non 
Homogenous Linear partial Differential 
Equations by Laplace transformations 


jULujVI 

olioiS AjjjIII 4ul£ / AjlaIa 




^ (JA *1% ^iC> ^ s j 


4 <ut il Vta^Ulj 4-uj^il AjIo^JI Aulualil!) 


iUuUaUll U^all <yfeiil 4^L«U*JI dfrUll J* 



Solving Homogenous and Non 
Homogenous Linear partial Differential 
Equations by Laplace transformations 

JfrluL » 1 \ jliulVI 


A«»rk<a (jiC- 

diLulI AjjjoII AalS / Ai j£J| AjLala. 


Cj|j ^ , „‘1 ri,„‘l-v~i^U AalaiJl AajjaJl AaLialalSI Ja. jA dtxgl lAA ^ U3AA 

QjdUi cIlJLjaj ^IAxLujIj AajliLSl cj^UU^i) 
(jiaauJ Ajuiali. AJjIxa!| ,jj£j Laiic , ^ 
.hjui ^U Vu ti) jjJj . Y 


Abstract 

Our aim in this paper is to solve the homogenous and 
nonhomogenous linear partial differential equations with 
constant coefficients by using Laplace transformations :- 

1. with some initial conditions. 

2. Without any conditions. 

Introduction 

Laplace Transformation [5] is an important 
transformations for solving linear ordinary and partial differential 
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equations with constant coefficients which subjected to some 
initial or boundary conditions. 

In [3] Mohammed, A. H was suggested different method 
from the usual method about solving linear ordinary differential 
equations by using Laplace transformations (L.T) with constant 
coefficients and submissive for initial conditions, also in [2] 
Mohammed, A. H suggested another method for solving linear 
ordinary differential equations with constant coefficients and 
without using any conditions by method of Laplace 
transformations also. In our work, we generalized the above two 
methods for solving the linear partial differential equations with 
constant coefficients . 

1. with some initial conditions 

2. without using any conditions. 

In [3] Mohammed, A. H supposed a linear ordinary 
differential equation of order (n) with constant coefficients and 
due to certain initial conditions, which general form can be 
written as follows 


a 0 y (n) + a 1 y (n L> + ... + a n y = f(x) 


, 0 - 1 ) 


Where f(x) is continuous function whose Laplace 
transformation can be determined and y(0),...,y ( " _1) (0) are 
defined, after he taking (L.T) to both sides to the above equation 
he put 


L(y) = 


h(s) 
k(s ) 


Where h and k are polynomials on (s), such that the degree 
of (h) is less than the degree of (k) and the polynomial (k) with 
known prime factors and by taking inverse Laplace 
transformation ( L~\T ) to the last equation, he put the general 
solution as follows: 


y = A 0 k 0 (x)+... + A m k m (x) 
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Such that k 0 ,...,k m are function of x ( m equal to the degree 
of k(s)) and that A 0 ,...,A m are constants. To find the values of 
A,. (i = 0 he substituted the initial conditions, one of them is 
y(0) and by taking derivatives of the last equation (m) times he 
got system of linear equations and by solving this linear system 
he obtained the values of A, (i = 0 and so by this method he 
obtained the solution of the required differential equation. 

In [2] Mohammed , A. H also, supposed the above 
(L.O.D.E) but without subjected to any initial conditions he took 
L.T to both sides, he found 


(a 0 s" + Oji '" -1 + ... + a n )h(s ) 

Where h(s) is polynomial of (s) represents denominator of 
L.T of the function f(x) and k(s) is also a polynomial of (s) and 
its degree smaller than the degree of the product of (a 0 s n + ... + a n ) 

and h(s) and after he took ( L~ l .T ) to the above equation he found 
that 

y = A Si (*) + - + A„g„ (x) + (x) + ... + B\ (x) 

Where A/s and B.'s are constants and g/s and h- A are 
functions of x. The number of constants B t 's and the number of 
h t 's are equal to the degree of h(s) which supposed to be ( r). He 
noted that the order of (O.D.E) is (n) , therefore its general 
solution must contain (n) constants only, but the last solution 
contains (n+r) constants and to solve this problem, he eliminated 
(Bj's) by obtaining whose values by substituting the last solution 

on (O.D.E), so he found a solution contains only (n) constants 
and by this method he got the general solution to the above 
(O.D.E) without using any initial conditions Mohammed, A. H 
[2] , [3] noted that the differ between his method and the usual 
method of (L.T) is that the usual method depends on the 
numerator, but his methods are independent on numerator. 
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Definition: The Laplace transform of a function f(x) is denoted 
by L{f(x)} and is defined as the integral of e fix) 
between the limits x=0 and x = °° , i.e 

L{fix)}=]e~ a fix)dix)=Fis) 

0 

The constant parameter (s) is assumed to be positive and 
large enough to ensure that the product e~ sx f(x ) converges to zero 
as x -> oo , for most common functions f(x) see [4]. 

Definition: Laplace transformation for a function u(x,t) is 
defined by the following form 

L{u(x,t)}= J e~ st u(x,t ) dt = v(x,s);s > 0 
0 

From above definition, we get the following 

L{u t } = sv(x, s ) - u(x, 0) 

L(u n } = s 1 2 v(x, s ) - su(x, 0) - u t (x,0) 

L{uJ = Dv(x,s);D = -^~ 

clx 

L{u xx } = D 2 v(x,s)-D 2 = ~~~2 

dx 

Note: For properties of Laplace transformations and it's 
inverse, see [1], [4]. 

The Suggested methods about partial Differential 
€quations:- 

1. The first method with some initial conditions, this method is 
summarized as follows :- 

Suppose we have a P.D.E of order (m+n) with constant 
coefficients due to certain initial conditions, which general form 
can be written as follows : 
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y r)/ 'c)x J 


= f(x,t ) ... (1) 


II 

i=0 7=0 

Where a. are constants, f(x,t) is a continuous function 


whose (L.T) can be determined, and some initial conditions 
u(x,0),u t (x,0),..., are defined. 


To find the solution of the P.D.E (1) we can take (L.T) to 
both sides of it and after substituting the conditions, we obtain an 
O.D.E of order (n) . 

The general form of it is 

a 0 v {n) +a 1 v < " _1) + ... + a n v = g(x,s) (2) 


Such that a 0 ,a 1 ,...,a n depend on (a. ) and (s) 

The later O.D.E its dependent variable is (v) and the 
independent variable is (x) , such that v=L {u} . 

After this, we write equation (2) by (D = — ) operator. 

dx 

Thus we get 

( a 0 D " + a,/)" -1 + ... + a n )v = g(x,s ) (3) 


From equation (3) , we get 


1 

N(D,s) 


■ g(x,s ) 


Where N(D,s) is a function of D and s. 
And after solving equation (4) , we get 

g(-r^) 

N(s) 


(4) 

(5) 


Where g(x,s) is a function of (x) and (s) and may be equal 
to g(x,s) or not, and N(s) is a function of (s). 

Now, since u = L~ l (v) , then 


u = L 1 { 


N(s) 


(6) 


Equation (6) represents the solution of the P.D.E (1), 
which we can be written by the form: 
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u = A 0 (x)k 0 (t) + A l (x)k 1 (t) + ...+A m (x)k m (t) (7) 

Such that k 0 ,k l ,...,k m are functions of (t) and that 
A 0 ,A l ,...,A m are functions of (x) , whose numbers equal to the 
degree of N(s) . 

The Heat Equation [6] 

The heat equation (one dimensional heat equation and 
without heating source) has the form 

u, = a 2 u „ 

Where u(x,t) represents the temperature at the point (x) 
and time (t), 

U l (x,t) represents the average of temperature change, 

U x (x, t) represents the slope of temperature curve, 

U xt (x, t) represents the convex or concave of temperature curve, 

k 

a constant which is called diffusibility and it's equal — , where 

at 

k represents heat conductivity coefficients and a represent 
density, and all of them are constants. 

We shall show the above method by the following 
examples 

Example (I) : To solve the irreducible P.D.E (Heat Flow 
Equation) 

U, = U xx ( 8 ) 


With only one initial condition u(x,0) = 3 sin 2nx , we take 
(L.T) to both sides of equation (8) , then 

d 2 v 


-u(x, 0) + sv 
d 2 v 


dx~ 


dx 9 

-sv = -3sin2;rx 


( 9 ) 


From equation (9) , we have 

v = — — (-3 sin 2nx) 

D~ -s 
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3 sin Ik x 

Thus, v = r— 

s+Ak 2 

Since, u = L _1 (v) 

So, u = 3e~ ^ sin 2k x (represents the solution of (8)). 


Example (S) : To solve the non homogenous reducible with 
homogenous terms P.D.E 

u x + 4 u t = -8t . . . (10) 

With initial condition u(x,0) =0 , we take (L.T) for both 

sides , we get 

dv . 8 f 

+ 4sv = — 5 - ... ( 11 ) 

dx s 2 

From equation (1 1) we get 
1 __ 8 _ 

V ~{D + 4s) s 2 

-2 

Since, u = L“‘ {v} 

So, u = - 1 2 (represents the solution of (10)) 


Example (3): To solve the reducible P.D.E (wave Equation) 
with homogenous terms 
u „ = ••• (12) 

With only two initial conditions u(x, 0) =sin;cc = -u t (x, 0) we 

take (L.T) to both sides of equation (12) , we get 

d 2 v 2 x • 

— - — s v = 1 - 5 ) sin ;tjc 
dx 2 

Thus 


v = 



(1 - s) sin kx 


^-1 

— ; -sin^x 

5 +K~ 


O \ 
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Since u = L 1 (v) => m = sin /dAcos sin ;r t ) 

TC 

Represents the solution of equation (12). 

2. Now, we will explain how we can use Laplace transformation 
to solve L.P.D.E without using any initial and boundary 
conditions. 

We will show this method through out second order linear 
partial differential equation with constant coefficients which has 
the general form 

Au xx + Bu xt +C u tl + ju x + Eu t + Fu = f(x,t ) (13) 

Which doesn't obey to any initial and boundary conditions 
and we suppose that the (L.T) of f(x,t) is known. 

To solve equation (13) we take (L.T) to both sides, we get 

(AD 2 + BsD + Cs 2 + JD + Es + F)v = ' h ^ (14) 

K (s) 

Equation (14) is an O.D.E of second order with constant 
coefficients A,Bs,J,Cs 2 +Es + F where k^s) is a polynomial of (s) 
represents the denominator of L.T of the function f(x,t), and 
h x {x,s) is a function of (x) and (s) and it isn't necessary to know 
the terms of \(x ,s) , we only denoted it by this symbol[2],[3]. 

To solve equation (14), we have 

v = ~ ~ . h t (x, s ) 

L (s)(AD~ + BsD + JD + Cs 2 +Es + F ) 

Now, we'll substitute any constant say (b) instead of D, 

(because the (L.T) is known for (e fl \sinax,cosax,sinhax,coshux) 

and its multiplications of these functions, but if the function f(x, 

t) consists polynomial in (x) we applied the fact 

— — = (1 + D + D 2 +...) , 

1-D 
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see [1], the constant b may be equal to a or ai(i = V 3 !)), 

that is 


h 2 (x, s) 

k, ( s)(Ab 2 + Bsb + Jb + Cs 2 + Es + F) 


_ h 2 (x,s ) 
k 2 (s) 

Where h 2 (x, s ) is a function of (x) and (s) and may be equal 
to \ (x, s) or not, k 2 ( 5 ) is a polynomial of (s) 

Since v = L{u} 


^u = L 1 } thus we h a ve 

*2 CO 

u = ^ B i(x)fi (0 (15) 


Number of i depends on the degree of k 2 (s) , 5, (x) are 
functions of (x) and flit) are functions of (t), (u) may contain one 
or more of the functions 5, (x) and f] (, t ) . 

To find the forms of the functions B : (x) , we take 
derivatives of equation (15) w.r.t (x) and (t) to get u x ,u xx ,u t ,u xt 
and u tt . 

We substitute all these derivatives and u in equation (13) 
and by equality coefficients of both sides we get the values of 
B j (x),\/i , and by this method we obtain the solution of the P.D.E 
(13). 


Example 4: - To solve the non- homogenous reducible P.D.E 

u xx~ u t, +U x +u,=e x+y ... (16) 

By using L.T and without using any conditions, we take 
L.T to both sides, we get 


(D 2 +D-s 2 +s)v = 


h{x, s) 
s - 1 


=> V = 


I 

(s-l)(D 2 +D-s 2 +s) 


h(x,s) 


(h consists e x ) 
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To find the solution of (v) we substitute (1) instead of the 
operator(D) (because the cofactor of a: for the function e x equal 
to one), so, we have 

_ h (x, s ) 

V ~ (s-l)(s-2)(s + l) 

_A(x) | B(x) | C(x) 

5-1 s-2 S + l 

Since u = L _1 {v} =^> u =A(x)e l +B(x)e 21 +C(x)e”‘ 

Where A, B and C are functions of x, and to find the forms 
of them, we must find u x ,u xx ,u t ,u lt and after finding them, we 
substitute them in P.D.E (16), thus we have 

A{x) = (c l +c 1 e ) + -^e x 

B(x) = c 3 e~ lx + c 4 e x 
C(x) = B(x) 

=> u = c^e’ + c 2 e~ x+t + ^ e x+l + {c 3 e~ lx + c 4 e x ){e 2t + e ~‘ ) 

Where ci,c 2 ,c 3 and c 4 are constants. 

Example 5 To solve the homogenous irreducible P.D.E 
u xx + u n + 2 u xt + 2 u x + 2 u t + u = 0 (17) 

We take (L.T) to both sides, we get 
(Z) + 2 sD + 2 D + s + 2s + 1) v = h(x, s) 

=> v = — — - ( ' X ’ ^ = — - — 5 - , (we substitute (a) 

a + s~ + 2s(ci + 1) + 2a + 1 (s + (a + 1))“ 

instead of operator D) 

A(x) | B(x) 

5 + (a + l) (.v + (a + l)) 2 
Since u = L~ l {v} 


0 t 
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=^> u = A<? _(fl+1) ' + 2?t <? _(fl+1) ' ; A and B are functions of (x). 

Now, to find the forms of the functions A and B , we must 
find u x , u xx , u t , u tt and u xt and after finding them, we substitute 
them in P.D.E (17), thus we find 

A(x) = (c t x + c 2 )e a ' - c 3 x 2 e ax 

B(x) = (c 3 x + c 4 )e ax ; C[ , c 2 , c 3 and c 4 are cons tan ts 
=> u = (c 1 x + c 2 - c 3 x 2 )e ax ~ (a+1)t + (c 3 x + c 4 )te ax ~ {a+1)t 
Where c 1 ,c 2 ,c 3 ,c 4 and a are constants. 
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